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Abstract 

We numerically calculate the first few eigenvalues of the perturbations 
of self-similar solutions of the spherically symmetric co-rotational SU(2) 
sigma-model on Minkowski space. 

1 Introduction 

In the recent past considerable interest focused on the study of nonlinear 
evolution equations that show "blow-up" for solutions in a finite time. 
Prominent candidates are self-similar solutions which were found to play 
an essential role in critical collapse phenomena. Critical solutions have 
the property that they lie at the boundary of two possible end states of 
dynamical evolutions: Typically, sufficiently large initial data will lead to 
the formation of some kind of singularity, while small data disperse. By 
fine-tuning the initial data to this threshold, the dynamics shows universal 
behaviour which is governed by a "critical solution". These solutions 
are characterized by the existence of a single unstable mode in linear 
perturbation analysis. While it is relatively easy to obtain the eigenvalue 
of the unstable mode by numerical analysis, to calculate the stable ones 
may be more subtle. However, the knowledge of the stable modes is 
important because they are responsible for the process of how fine-tuned 
data are attracted toward the critical solution. 

In this note we study equivariant maps of the SU(2) sigma-model (wave 
maps) from 3+1 dimensional Minkowski spacetime into the three-sphere. 
This system admits a discrete sequence of self-similar solutions and has 
been studied as a toy model for critical collapse behaviour [5J. Linear 
perturbation analysis shows that the "ground state" is stable, while the 
higher excitations are unstable. The first excitation acts as the critical 
solution between collapse to the ground state and dispersion. Bizori jl] 
has given a beautiful algorithm how to obtain analytically all eigenvalues 
for the linear perturbations of the ground state. However, for the higher 
excitations, because they are known only numerically, one has to rely on 
numerical tools. 
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We apply two different numerical methods to calculate the first few 
eigenvalues of the ground state and the first excitation i.e. the critical 
solution. One method is by shooting and matching an the other by time 
evolution. Our analysis shows that the "shooting" method, frequently 
applied for finding the largest eigenvalues may fail for the smaller values. 

For the ground state our results confirm the analytical results. This 
is important because, as pointed out by Bizoh, the eigenvalue problem 
is not standard. Our results show that the analyticity condition for the 
singular boundary value problem chosen by Bizoh are in accordance with 
the time evolution of the system. In addition we find for the critical 
solution, besides of the unstable and the gauge mode, the first two stable 
eigenvalues. 



2 The SU(2) a-model 

2.1 Definition of the model 

A smooth mapping U : (M,g) — > (S 3 , G) from a spacetime M with metric 
g into the three-sphere is called SU(2) cr-model if it is a critical point of 
the action 

f g ab d a U A d b U B G A B. (1) 

J M 

Let M be Minkowski spacetime. We use polar coordinates (t, r, 8, ip) on 
M and standard coordinates (u, ©,<!>) on S 3 . In these coordinates the 
metrics are given by 

g = ~dt 2 + dr 2 + r 2 (d6 2 + sm 2 {9)d(p 2 ) 

G = du + sin 2 (u)(dQ 2 + sin 2 (e)d$ 2 ) (2) 
We restrict ourselves to maps of the form 

U(t,r,$, V ) = {u(t,r),0,<p) (3) 

which are called spherically symmetric and co-rotational. Under these 
assumptions the Euler-Lagrange equations associated with the action Q 
reduce to the single semilinear wave equation 

2 sin(2u) . , 

Utt - Urr - -U r H ^ = 0. (4) 

2.2 Self-similar solutions 

A self-similar solution of eq. |gj is a solution of the form 

"(*'»•) = /(T3t) W 

for some constant T > 0. We substitute the ansatz ® into eq. IjlJ and 
obtain 

r + V-^Vo ( 6 ) 
p p 2 (i-p 2 ) 

where p := jf—^ and ' := Note that the domain p G [0, 1] is the back- 
ward lightcone of the point (T, 0) which is called the singularity because 
u r (t, 0) grows unbounded for t — > T-. Therefore T is called the blow-up 
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time. Eq. (|SJ has two singular points (p — 0, 1) and therefore we have to 
require the boundary conditions 

/(0) = and /(l) = | (7) 

to ensure smoothness of solutions. Bizoh has shown existence of a 
countable family {fo, ,fi, . . .} of smooth solutions of eq. © which satisfy 
the boundary conditions Q. The so-called ground state fo is known in 
closed form and given by 

fo (p) — 2 arctan p. (8) 

The higher excitations /i, /a, . . . can be constructed numerically using a 
standard shooting technique (cf. 



2.3 Perturbations 

We are interested in linear perturbations around the solutions fo and /i 
which play a crucial role in dynamical evolution. We introduce so-called 
adapted coordinates (r, p) defined by 

T = -log(T-i)andp=^-^. (9) 

These coordinates cover the domain t < T (the singularity is shifted to 
r — > oo). In the new coordinates eq. Q becomes 

u TT + 2pu Tp - (1 - p 2 )u pp + u T - 2 1 - P u p + 8m @ M ) = o. (10) 

p p 2 

It is important to note that the class of self-similar solutions which blow 
up as t — > T— are static (i.e. independent of r) in these coordinates. 
Inserting 

u(r,p) = fn(p) +w(r,p) (11) 

into eq. (1 M ) I and neglecting terms of order w 2 we obtain a time evolution 
equation for linear perturbations around the n-th self-similar solution f n . 

w TT + 2pw Tp - (1 - p' 2 )w P p + w T - 2- —w p + 2cos ( 2 /") w — Q (12) 

p p 2 

To ensure regularity we have to require the boundary conditions 

w(t, 0)=0 and w pp (t, 0) = 0. (13) 

We are looking for perturbations of the form 

w(t, p) = e Xr v\(p) (14) 

where A may be complex in general. A is said to be an eigenvalue and 
v\ an eigenmode or eigenfunction. Like in the finite-dimensional case one 
expects that these eigenmodes determine the dynamics in the neighbour- 
hood of the solution /„ in phase space. 

To obtain an ordinary differential equation for the eigenmodes v\ we 
insert ansatz l|14|l into eq. I|12|l . 

„, (A + l)p 2 -l , A(A + l)p 2 +2cos(2/4p)) 

MP) + %(P + D(P-1) ^ (P) + P 2 (P + 1)(P-1) MP) = °' 

(15) 
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This equation has again two singular points. The regularity conditions 
are given by 

«a(0) = 0, wi(0) = o 
VA (l) = l,^(l) = -^±^ (16) 

where a £ R is a free parameter and we have normalized the eigenfunctions 
to have V\(l) = 1. 



3 Numerical calculation of eigenvalues 
3.1 Calculating eigenvalues by time evolution 

We consider the Cauchy problem eq. 1121 together with initial data 
w(0,p) = /(p) and w T {0,p) = g{p) inside the backward lightcone of the 
singularity, that is for p £ [0, 1]. By setting u\ = w, U2 — w T , 113 — w p 
and u — (iti, 1*2,1*3) this problem can be written as a first order system 

-^-u = Lu (17) 

(IT 

with initial conditions lti(0,p) = f(p), 1*2(0, p) = g(p) and i*3(0,p) = 
/'(p). L is the linear differential operator following from eq. il'It . If the 
Cauchy problem is well-posed one can introduce a one-parameter family 
of bounded linear operators S(t) for r > such that the solution to initial 
data it(0, •) is given by 

u(r, •) = S(t>(0, •)■ (18) 

We are looking for special solutions defined by eq. 11411 . Inserting this 
ansatz into eq. I|18^ one obtains 

S(t)(v\, \v x ,v' x ) = e XT (v\,\v x ,v' x ) (19) 

for an eigenmode v\ with eigenvalue A. For fixed r we define an inner 
product by 



rl 3 

(u(t, -)\v(t, ■)) := / ^Uj{T,p)vj{T,p)dp 
Jo j=1 



(20) 



and a norm by 



Hr, .)|| :=V<«(t,.)Kt,.)>. (21) 
If v\ is an eigenmode with eigenvalue A then we have 

\\S(r)(v x , A«a,«a)|| = |e Ar 1 1| K , \v x , v' x )\\ = Ce RcAT (22) 

where C is constant. Now we make the following assumptions: 

• The eigenvalues {Ao, Ai, . . .} form a discrete set 

• Arbitrary initial data can be decomposed into a sum of eigenmodes 

Then we can write (real- valued) initial data u(0, •) as 

00 00 
«(0,p) = ^(c k ipk(p) +c k ip k (p)) = 2^Re(c fe </?fe(p)) (23) 

k=0 fe=0 

where ip k := (v\ k ,AkV\ k ,v' Xk ) and v\ k is an eigenmode with eigenvalue 
A fc . 
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For simplicity we write fik '■= ReA^, uj^ := ImAfc and without loss of 
generality we assume that y,o > fix > , . .. Applying the time evolution 
operator S(t) to initial data ti(0, •) yields 

oo oo 

S(t)u(0, •) = 2^Re( Cfe e A *>fc) = 2^ e ^ T Re(c fc e l ^> fc ). (24) 

k=0 k-0 

So for large r and > Mi we have 

\\S(r)u(0,-)\\^e^g(r) (25) 

where g(r) := ||Re(coe*"° T <po)||- If there are more eigenvalues with the 
same real part po then the function g in 125H is more complicated but the 
result is essentially the same. 

For computing the real part no of the dominant eigenvalue one simply 
evolves arbitrary initial data and calculates the norm 12 U . As long as 
the expansion coefficient Co in i'2'M does not vanish, the logarithm of the 
norm will behave as hot + log(g(r)) for large r and by least-square fitting 
one can read off /io- Note that g is constant if and only if Ao is real. 

To calculate real parts of other eigenvalues we use the orthogonaliza- 
tion method described in 0. Here we only sketch the basic idea. Take lin- 
early independent initial data <f> and ip. Let P^ denote the orthogonal pro- 
jection on {4>} ± ■ One can readily show that Ps(T)tf>S(T)P,t, = -Ps( T )^5(r). 
Now we define a third function \ '■— ip + a 4> where a is some constant. 
Clearly we have P$x = P$ip. Suppose we choose the constant a in such a 
way that \ nas the expansion 

oo 

X = 2^Re(c fc¥ 3 fc ). (26) 

k = l 

Note that the sum starts at k = 1! For simplicity we assume that Ao and 
Ai are real. Then, if Ao > Ai > ReA2, 

S(t)x = 2e AlT Re(ci^i) + . . . (27) 

Unfortunately we do not know the right a which is required to obtain 
expansion (12(H . But surprisingly this is not necessary since 

P s(t)4 ,S(t)<P = P s(r)r „S(T)P^ = P s{t)4> S{t)P4, X = 

Ps( T)< t>S{T)x = e AlT P s(rW Re(c m ) + . . . (28) 
Observe that for real Ao we have 

Ps{t)4, ~ P-Rc(c ovo ) (29) 

for large t. It follows that 

||P SCT) ^(r)V|| ~ e Air (30) 

for large r and we can read off Ai. In practical application the orthog- 
onal projection is performed at every timestep (see |3] for details). The 
probably most important feature of this method is that it works with 
an arbitrary inner product. The eigenfunctions ipk do not have to be 
orthogonal. 
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3.2 Calculating eigenvalues by shooting 

Another possiblity for calculating eigenvalues is to use eq. (1151) . One has 
to solve an ODE boundary value problem where the boundary conditions 
are given by the regularity requirements . This problem can be solved 
using a standard shooting to a fitting point technique. Although this 
method is very accurate it has at least two serious disadvantages: 

• It is possible to "miss" some eigenvalues because one has to start 
with sufficiently good initial guesses of the shooting parameters. 

• The method cannot be expected to work for large negative eigenval- 
ues (see 

Nevertheless it is a good check for results obtained by time evolution. 

4 Numerical results 

We have calculated four eigenvalues of the ground state and the first exci- 
tation using the orthogonalization method. The first negative eigenvalue 
can be obtained by shooting as well. The positive eigenvalues have been 
computed before in via shooting. Bizon A. has also calculated a large 
number of negative eigenvalues for the ground state using a completely 
different method. His results are in good agreement with ours. Further- 
more, in [I] it has been shown analytically that A = —2 is an eigenvalue 
of the ground state. 

The eigenvalue A = 1 corresponds to a so-called gauge mode which 
is associated with the freedom of choosing the blow-up time T in the 
coordinate transformation © (see [5] for a more detailed explanation). 
This gauge mode has no physical relevance. 

We summarize our results in the following two tables. 



Eigenmode 


Time Evolution 


Shooting 


Gauge 


1 


1 


First stable 


-0.5424 


-0.54246 


Second stable 


-2.00 




Third stable 


-3.3 




Table 1: Eigenvalues of the ground state 


Eigenmode 


Time Evolution 


Shooting 


Unstable 


6.3336 


6.333625 


Gauge 


1 


1 


First stable 


-0.518 


-0.5186 


Second stable 


-1.7 





Table 2: Eigenvalues of the first excitation 

Remarkably we do not see any oscillations and therefore all eigenvalues 
seem to be real. 
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A Numerical integration of the evolution 
equation 

A.l Integration scheme 

For integrating eq. (1171 we apply the second-order characteristic method 
used e.g. in pj- We write eq. I|17|l as 

du A du . 
— =A— + Bu 31 
or dp 

where A and B are 3x3 matrices depending on p. The canonical decom- 
position of the matrix A is given by 

A = X-A- +X+A+ (32) 

where A_, A+ are the two non-zero eigenvalues of A. A- and A+ are the 
orthogonal projections on the eigenspaces. The (signs of the) eigenvalues 
of A define the directions of the characteristics of the system (1311) . We 
discretize the system (1311 in space by 

(fit 

^- = \-A-D{\-)u + \+A+D{\ + )u + Bu (33) 

where D(\) is a finite difference operator producing either 3-point left- 
sided derivatives or 3-point right-sided derivatives depending on the sign 
of A. For the discretization in time we use a Runge-Kutta-like scheme (cf. 
0). 



A. 2 Boundary conditions 

At the boundary p — 1 we do not have to impose any boundary condition 
because it turns out that A± < for p > 1. Therefore we do not need 
right-sided derivatives here. Physically this is clear since p = 1 is the 
past lightcone of the singularity and therefore no information can come 
in from outside. At p = 1 — Ap (where Ap denotes the mesh of our 
grid) we do need a right-sided derivative and since we are using 3-point 
approximations we need one point outside the lightcone p = 1. Therefore 
we integrate the equation on p £ [0, 1 + Ap]. 

At the center p = we have the regularity conditions (cf. 1131 ) 

Ui (r, 0)s0 and ^(t,0)sO. (34) 

To impose these boundary conditions we define ghost points by 

M i( r .-p) = - M i( r .p) (35) 
and use the same integration scheme as in the interior. 
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